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Transformation de Fourier discrète et algorithmique. Théorème de Courant-Fischer et applications. 

Corrigé par M.TARQI 1 http ://alkendy.xl0.mx 

Premier problème 

l ere partie 

Quelques propriétés de ( 1> P 

1.1 


1.1.1 <3? p (x) = 0 si et seulement si y^ = 0, pour tout k = 0,1 — 1, donc 3> p (x) = 0 si 

et seulement si P x annule les w p , donc P x est un polynôme nulle (deg P x < p — 1 et les 
scalaires w p sont deux à deux distincts ). 

1.1.2 II est clair que <L p est un endomorphisme de C p , de plus si <L p (x) = 0 si et seulement si 
P x = 0 ou encore (xq, x\ . x p _i) = x = 0, donc <L p est un automorphisme de C p . 

1.2 


1.2.1 Si x 

rriij 


Ci, alors P x = X 1 et donc <& p (ej) = ( 1, , w p . . . , w p , 1 ' 1 *) et par conséquent 

(r^p )o<i.]<p 1 ■ Ainsi 


( 1 1 

1 w p 


1 \ 


V 1 


W' 


( P-1 


W. 


(p- 1) 2 




1.2.2 La matrice M est une matrice de Vandermonde, et comme les nombres 1, iu p . w p . vfP 1 
sont deux à deux distincts, alors M est inversible, donc on retrouve le fait que <1> P est un 
automorphisme de CL 
1.3 


1.3.1 Les (w^ )/ ;Æ m étant les racines p-ièmes de l'unité, donc en tenant compte de la relation entre 

p- 1 ., fc p- 1 

les coefficients et les racines, w p J = 0 si p ne divise pas i — j et w p ' ' = P si p 


k = 0 


k = 0 


P - 1 


divise i — j. Ainsi si (i, j) sont dans {0, 1, 2, ...,p — 1}, w p = pàij> °ù Sij désigne le 


k = 0 


symbole de Kronecker. 

1.3.2 Soit (xo,xi, x p _i) G C p et (yo, y ±, ..., y p _ i) = <& p (x), alors on peut écrire : 

p— 1 p— 1 p— 1 p— 1 p— 1 \ 


X j W p 


J^\yk\ 2 = p x( w p)p*( w p) = ,jhw p' 5Z 

/ c =0 / c =0 / c =0 2—0 2—0 

p— 1 / p— 1 \ p- 1 p— 1 

= S 

k= 0 \i=0 i^j 

p— 1 p— 1 




I - 


p- 1 




w. 


k{i-j) 


k=0 2—0 




/c=0 


pX^i 2 + 0 


i=0 


i=0 


1 Si vous avez des critiques ou des encouragements à formuler sur le contenu, n'hésitez pas à nous en faire part, et 
surtout n'hésitez pas de me signaler les erreurs rencontrées. 
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1.4 


1.4.1 Notons MM = (cij)o<ij< p - 1 , alors pour tout ( i,j ), on a : 


p — 1 p — 1 


c *i = 2^ w p w p 

k = 0 /c=0 


p, SI Z = J 

0, sinon 


D'où: 


MM = p/p. 


1.4.2 La relation MM = pl p montre aussi que M est inversible et que M 1 = -M. 


Remarque : Si on pose M = M Wrt , alors M 1 = - AL _i. 

x A ” p Wp 


2 ème partie 

Un peu d'algorithmique 


2.1 


P 

2.1.1 D'après l'algorithme, on a pour tout & = 0, 1, — 1 : 

F k = w pi k 
a-k = Pk + Wplk 
Ac+f = Afc - ^p7fc 

2.1.2 Par définition de <f? a , on a pour tout k = 0,1, ...,p — 1 : 


p - 1 


f- 1 


f- 1 


a k = Y2 a i W p j = 22 a ^ W p 2j + 22 a ‘2j+\W 


k(2j+l) 


1=0 


3=0 


3=0 


E- 1 
2 X 


£-1 
2 X 


Ceci s'écrit aussi ccfc = a 2 jW^ + w p ^ a, 2 j+\w k J . 


Finalement si 0 < k < 


3=0 

P 


3=0 


et 


et donc 


«fc = Pk + Wp7fc 

E ( fc +§)1 k 

~ W p2^ °2i+! 
1=0 2 1=0 


2-1 
2 1 


2-1 
2 A 


(*+§)! 
tc E 
2 


Ofc+E = Afc - Wp7fc- 


Donc l'algorithme de la question 2.1 permet de calculer de proche en proche toutes les 
composantes de 4> p (a). 


2.2 
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2.2.1 Dans cette question on supposera que les différentes puissances de w p ont été calculées. 
On a «i) = (ao+ai, ao — ai), donc si = 2 et r\ = 0. La calcul de 4>2«(a) fait intervenir 

deux transformations de Fourier de taille | à savoir ( l>v (b) et ( \> v (c), ainsi | multiplications 
complexes et p additions complexes. On a donc les formules de récurrences suivantes : 

j s n = 2s n -\ + p = 2s n -i + 2 n 
| r n = 2r n _i + | = 2r n _i + 2 n_1 , n > 1. 

2.2.1 Notons u n = s n + r n le nombre total des additions et de multiplications complexes néces- 
saires au calcul de <ï> 2 n (cl), alors on a : 

'«i =2 et u n = 2u n -i + 3 x 2 n ~ 1 , n > 1. 

Ainsi u n = 2 n + 3 (n - 1)2”“ 1 =p + ^p(log 2 p - 1). 

2.3 

2.3.1 Pour évaluer le polynôme en A on pose : 

Up— 1 — dp— 1 > Up — 2 — o,p — 2 -F A yp—\,...,yp—] i — dp— k + Ay n _fc_|_ pour k — 3, ■■•,p 1. 

Alors yo = P a ( A) et donc on a besoin de 2 (p — 1) opérations : p — 1 additions et p — 1 

multiplications. 

2.3.1 Notons v v le nombre total des additions et de multiplications complexes nécessaires au 
calcul de $ p (o) par la méthode de Horner. Pour calculer (a) on aura besoin de p valeurs 
Vk = Pa(w k p ), k = 0, 1, ...,p — 1, donc le calcul de <!> 2 n (d) demande 2 p(p — 1) opérations, soit 
v p = 2 p 2 — 2p. 

2.4 La méthode de division décrite dans la question 2.1.2 demande u n opérations avec 

3 

Un =P + ~p(log 2 p- 1) = 0(p\0g 2 p), 

alors que la méthode d'Hôrner demande v n opérations avec 

v p = 2 p 2 -2 p = 0(p 2 ). 

Conclusion : le calcul par la première méthode est donc très nettement plus rapide que celui 
par la deuxième méthode. 


Deuxième problème 

l ere partie 

Théorème de Courant-Fischer 

1.1 Puisque A est symétrique, l'endomorphisme f,\ est symétrique. D'après le théorème spectrale, 
toute matrice symétrique réelle est orthogonalement diagonalisable, autrement dit M" admet 
une base orthonormée de vecteurs propres de / a ■ 

1.2 


1.2.1 Pour tout fc = 1,2, ...,n, on a /^(efc ) = A^e^, et donc RA(ek) 
Afc- 


(/A(efc)|e fc ) 

(d-k \dk) 


(^kd-k \d-k) 
(dk | dk) 
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k 

1.2.2 Soit v = ViGi € Vk un vecteur non nul, alors : 

i = 1 


(Ia(v) | u ) = y XiVi \ Viei ) = XiV t - Xk Yl v i 


i— 1 i— 1 


2—1 


2—1 


et donc 


Ra(v) = tMM < v 

(v\v) 


D'autre part, on sait que R A (e k ) = X k avec e k G 14\{Q}, donc A/, = max R A (u ) 

W6l4\{0} 


1.3 


1.3.1 Si dim(F n Vect (e^ . ..., e n )) = 0, alors on aura : 

dim ( F + Vect(efc, ..., e n ) = dirn F + dim Vect(efc, ..., e n )) = k + n — k + 1 = n + 1, 
et ceci est absurde, donc nécessairement dim(F fl Vect (e k , ..., e n ) ) > 1. 

n 

1.3.2 Soit w=J2 w i e i G F fl Vect(efc, ..., e n ) non nul, alors alors : 

i=k 

n n \ n n 

(f A (w)\w) = Y X i W i I Wiei ) = X] XiW i ~ Xk 

i=k i—k / i—k i=k 

et donc 

r, , A Ua(w)\w) ^ s 
Ra[w) = — , i A — > Xk- 
(w\w) 

1.4 D'après ce qui précède, pour tout F G Tk, il existe te G F non nul tel que Fa (te) > Xk, donc 

pour tout F G Ffc max R A ( te) > Xk- Or il existe F G Tk, F = 14 tel que Afc = max Ra(v), 
w&f\{o } ^ ueUfcVo} 

on en déduit l'égalité demandée : 


Xk 


min 

F&F k 


max R a (v ) 
»6F\{0} 


2 ème partie 

Continuité de dérivabilité des valeurs propres d'une applications matricielle 


2.1 


Par définition de la norme subordonnée et l'inégalité de Cauchy-Schwarz, on a, pour tout v G 

■XAn,i . 

|(Ce|e)| < ||Ce||||e|| 2 < ||C|||M|| 


et donc si v / 0, 


(Cv\v) 

(e|e) 


< HCl 


2.2 


2.2.1 Soit v G Vk(t 0 ) non nul, alors on a : 

RA(t)(v) = RA(t)-A(to)( v ) + R A(to)(v) < Afc(to) + sup{R A (t)-A(t 0 )(v)/v € V k (t 0 )\{0}} 
En particulier, si v G V k (to) fl Vect(efc(t), ...,e n (t)) on aura A k(t) < RA(t)( v ) et donc 
A k(t) < X k (t 0 ) + sup{R A (t)-A(t 0 )(v)/v G 14(t 0 )\{0}}. 
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2.3 


2.2.2 Soit t, ta G I. Puisque t et tç, jouent le même rôle, l'inégalité précédente de la question peut 
s'écrire aussi : 

|A k(t) ~ A fc (t 0 )| < |sup {R A (t)-A(t 0 )(v)/v G Vfc(t 0 )\{0}}| 

mais la question 2.1 de cette partie montre que : 

\RA(t)-A(to)(v)\ < 11-4(0 — ^.(^o)ll- 

D'où 

V(t, to) G 1 2 , |A k(t) - A fc (t 0 )| < p(0 - Al(f 0 )||, 

et la continuité de A entraîne celle de A 

2.3.1 L'application A est C 1 sur R si et seulement si les applications a et b sont de classe C l sur 
R. 

• a est de classe C 1 sur R* et lim ait) = 0, donc a est continue sur R ; 

o 

• Vf / 0, a\t) = ^ e & cos (^j + sin ; 

• lim — — = 0 et lim a'{t) = lim x 3 e _æ2 (2xcosx + sinx) = 0. 

t— >0 t t — >0 |x|— >oo 

Ainsi a est de classe C 1 sur R, le même raisonnement se fait pour montrer que l'application 
b est de classe C 1 sur R. 

2.3.1 Pour tout t G R, A(t) est symétrique donc est daigonalisable sur R, Notons Ai (t) et A2 (t) 
les valeurs propres de A(t) et v\ (t). V 2 U) deux valeurs propres respectivement associées. 
Si t = 0 A(t) = 0 et Ai(0) = A 2 (0) = 0. 

Soit maintenant t / 0, le polynôme caractéristique de A(t) est 

A 2 -a 2 lt) -b 2 lt) =X 2 -e^, 

-1 -1 

et donc les valeurs propres sont Ai (t) = e l 2 " et A2(f) = — e"? 3 ". 

• Le sous espace associé à Ai (t) est engendré par l'application v\ (t) = 



si 


. 1 fl 
M — et ci 


2&7T \2kir 

Le sous espace associé à A2 (t) est engendré par l'application V2 (t) = 


/-«(IL 

l 2t. 

1 cos — . 

V V2 tj ) 


SI 


t + 


1 


et ci 


1 


(2 k + 1 ) 7r y (2 k + 1)tt / 

2.3.3 Supposons qu'une telle fonction e existe, posons e(t) = (ei (t), 62 (t)), alors on aura : 

Vf G R, A(t)e(t) = Ai (t)e(t) ou A(t)e(t) = A2 (t)e(t). 

SiA(t)e(t) = Ai (t)e(t), alors 

(*) Vf G R*, ei(t) = cos \ ei (t) + sin (- j e2(f). 
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Si t = 


1 


(2k + 1)tï 


, la relation (*) donne e\ 


1 


(2k + 1)7 r 

(1) ei (o) = o. 


Si t = 


(1 + 4/c)7T 


, on obtient e\ 


(1 + 4A:)7t 


= e 2 


= 0 et par continuité on aura : 


et de même raison on aura : 


(1 + 4£:)7t 


(2) ei(0)=e 2 (0). 


On déduit de (1) et (2) que ei(0) = e 2 (0) = 0 et ceci est absurde puisque e est un vecteur 
propre. 

La même conclusion dans le cas où A(t)e(t ) = À 2 (t)e(t). 

2.4 Soit t G I. M(t) étant symétrique réelle, donc elle est diagonalisable, soient Ai (t) et À 2 (t) ses 
valeurs propres 

a(t) b(t) A a c gon |. des applications de classe C 1 sur 
b(t) c(t) J 

I dans R. Le polynôme caractéristique de M (t) s'écrit donc : 


Posons, pour tout tG/, M (t) = 


XM(t)(X) = x 2 - (a(t) + c(t))X + (a(t)c(t) - b(t ) 2 ) 

Il admet pour discriminant 

A (t) = ( a(t ) + c(t)) 2 - 4 (a(t)c(t) - b(t ) 2 ) = ( a(t ) - c(t )) 2 + (2 b(t)) 2 

et d'après les donnés, il existe une application p : I — > R de classe C 1 telle que 

Vt € I, ( a(t ) - c(t)) 2 + (2b(t)) 2 = fi 2 (t), 

donc les valeurs de M(t ) sont Ai(t) = p(t) et A 2 (t) = — q(t), elles sont bien de classe C 1 sur I. 


2.5 


2.5.1 Si t = 0, B(t) = 0. Soit maintenant t / 0, alors le polynôme caractéristique de B(t) s'écrit : 

* 2 - *“ (2 + si “ (7) ) = (x - i 2 (2 + » (i) ) ) (x + ( 2 + sin Q) ) ) 

Donc le polynôme caractéristique de B(t) est scindé à racines simples, donc B(t) est dai- 
gonalisable dans A4 2 (M). 

2.5.2 B est de classe C 1 sur M si et seulement si b est C 1 sur M. Il est clair que b est continue sur M 

et dérivable sur R*, en 0, lim — — = 0 et lim b'(t) = 0. Donc B est C 1 sur M. 

t- 0 t t- 0 w 

Mais si on désigne par A (t) = t 2 ^2 + sin si t 7^ 0 et A(0) = 0, alors A est continue 

sur M, dérivable sur M* mais en 0, lim A '(t) n'existe pas puisque 


pour tout t non nul. 


A'(f) 


4 1 + 2 1 sin 




Fin de l'épreuve 
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